Search for: All records

Abstract Pulsar timing arrays (PTAs) are approaching the sensitivity required to resolve gravitational waves (GWs) from individual supermassive black hole (SMBH) binaries. However, the large uncertainty in source localization will make the identification of its host environment challenging. We show how to convert the posterior probability function of binary parameters inferred by GW analyses into distributions of apparent magnitudes of the host galaxy. We do so for a scenario in which the host environment is a regular early-type galaxy, and one in which it is an active galactic nucleus. We estimate the reach of PTAs in the near and intermediate future, and estimate whether the binary hosts will be detectable in all-sky electromagnetic (EM) surveys. A PTA with a baseline of 20 yr and 116 pulsars, resembling the upcoming data release of the International Pulsar Timing Array, can detect binaries out to a luminosity distance of 2 Gpc (corresponding to a redshift ofz ∼ 0.36), while a PTA with a baseline of 30 yr and 200 pulsars can reach out to distances slightly greater than 3 Gpc (z ∼ 0.53). We find that the host galaxies of all binaries detectable with a baseline of 20 yr are expected to be present in the Wide-field Infrared Survey Explorer and SuperCOSMOS surveys, if they lie outside the plane of the Milky Way. The Two Micron All Sky Survey becomes incomplete for hosts of binaries more massive than 10^9.8M_⊙at a luminosity distance greater than 1 Gpc. The EM surveys become slightly more incomplete when PTAs with longer baselines and therefore improved sensitivities are considered. 

Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes (SMBHs). These searches, across existing time-domain surveys, have produced several hundred candidates. The general stochastic variability of quasars, however, can masquerade as a false-positive periodic signal, especially when monitoring cadence and duration are limited. In this work, we predict the detectability of binary SMBHs in the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We apply computationally inexpensive sinusoidal curve fits to millions of simulated LSST Deep Drilling Field light curves of both single, isolated quasars and binary quasars. The period and phase of simulated binary signals can generally be disentangled from quasar variability. Binary amplitude is overestimated and poorly recovered for two-thirds of potential binaries due to quasar accretion variability. Quasars with strong intrinsic variability can obscure a binary signal too much for recovery. We also find that the most luminous quasars mimic current binary candidate light curves and their properties: The false-positive rates are 60% for these quasars. The reliable recovery of binary period and phase for a wide range of input binary LSST light curves is promising for multi-messenger characterization of binary SMBHs. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 mag will result in samples that are overwhelmed by false positives. This paper represents an important and computationally inexpensive way forward for understanding the true and false-positive rates for binary candidates identified by Rubin. 

Extremely variable quasars can also show strong changes in broad-line emission strength and are known as changing-look quasars (CLQs). To study the CLQ transition mechanism, we present a pilot sample of CLQs with X-ray observations in both the bright and faint states. From a sample of quasars with bright-state archival SDSS spectra and (Chandra or XMM-Newton) X-ray data, we identified five new CLQs via optical spectroscopic follow-up and then obtained new target-of-opportunity X-ray observations with Chandra. No strong absorption is detected in either the bright- or the faint-state X-ray spectra. The intrinsic X-ray flux generally changes along with the optical variability, and the X-ray power-law slope becomes harder in the faint state. Large-amplitude mid-infrared variability is detected in all five CLQs, and it echoes the variability in the optical with a time lag expected from the light-crossing time of the dusty torus for CLQs with robust lag measurements. The changing-obscuration model is not consistent with the observed X-ray spectra and spectral energy distribution changes seen in these CLQs. It is highly likely that the observed changes are due to the changing accretion rate of the supermassive black hole, so the multiwavelength emission varies accordingly, with promising analogies to the accretion states of X-ray binaries. 

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) atz= 0.359 with 153 spectroscopic epochs of data representing a 10 yr baseline (2013–2023). We split the baseline into two regimes based on the 3× flux increase in the light curve: a “low state” phase during the years 2013–2019 and a “high state” phase during the years 2022–2023. The velocity-resolved lag profiles (VRLPs) indicate that gas with different kinematics dominates the line emission in different states. The HβVRLP begins with a signature of inflow onto the broad-line region (BLR) in the low state, while in the high state it is flatter with less signature of inflow. The HαVRLP begins consistent with a virialized BLR in the low state, while in the high state shows a signature of inflow. The differences in the kinematics between the Balmer lines and between the low state and the high state suggests complex BLR dynamics. We find that the BLR radius and velocity (both FWHM andσ) do not obey a constant virial product throughout the monitoring period. We find that the BLR lags and continuum luminosity are correlated, consistent with rapid response of the BLR gas to the illuminating continuum. The BLR kinematic profile changes in unpredictable ways that are not related to continuum changes and reverberation lag. Our observations indicate that nonvirial kinematics can significantly contribute to observed line profiles, suggesting caution for black hole mass estimation in luminous and highly varying quasars like RM160. 

ABSTRACT Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g. Active Galactic Nuclei (AGNs) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >108M⊙) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multimessenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multimessenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly. 

Abstract Quasars atz≳ 1 most often have redshifts measured from rest-frame ultraviolet emission lines. One of the most common such lines, Civλ1549, shows blueshifts up to ≈5000 km s⁻¹and in rare cases even higher. This blueshifting results in highly uncertain redshifts when compared to redshift determinations from rest-frame optical emission lines, e.g., from the narrow [Oiii]λ5007 feature. We present spectroscopic measurements for 260 sources at 1.55 ≲z≲ 3.50 having −28.0 ≲M_i≲ − 30.0 mag from the Gemini Near Infrared Spectrograph–Distant Quasar Survey (GNIRS-DQS) catalog, augmenting the previous iteration, which contained 226 of the 260 sources whose measurements are improved upon in this work. We obtain reliable systemic redshifts based on [Oiii]λ5007 for a subset of 121 sources, which we use to calibrate prescriptions for correcting UV-based redshifts. These prescriptions are based on a regression analysis involving Civfull-width-at-half-maximum intensity and equivalent width, along with the UV continuum luminosity at a rest-frame wavelength of 1350 Å. Applying these corrections can improve the accuracy and the precision in the Civ-based redshift by up to ∼850 km s⁻¹and ∼150 km s⁻¹, respectively, which correspond to ∼8.5 and ∼1.5 Mpc in comoving distance atz= 2.5. Our prescriptions also improve the accuracy of the best available multifeature redshift determination algorithm by ∼100 km s⁻¹, indicating that the spectroscopic properties of the Civemission line can provide robust redshift estimates for high-redshift quasars. We discuss the prospects of our prescriptions for cosmological and quasar studies utilizing upcoming large spectroscopic surveys. 

Abstract “Changing-look” active galactic nuclei (CL-AGNs) challenge our basic ideas about the physics of accretion flows and circumnuclear gas around supermassive black holes. Using first-year Sloan Digital Sky Survey V (SDSS-V) repeated spectroscopy of nearly 29,000 previously known active galactic nuclei (AGNs), combined with dedicated follow-up spectroscopy, and publicly available optical light curves, we have identified 116 CL-AGNs where (at least) one broad emission line has essentially (dis-)appeared, as well as 88 other extremely variable systems. Our CL-AGN sample, with 107 newly identified cases, is the largest reported to date, and includes ∼0.4% of the AGNs reobserved in first-year SDSS-V operations. Among our CL-AGNs, 67% exhibit dimming while 33% exhibit brightening. Our sample probes extreme AGN spectral variability on months to decades timescales, including some cases of recurring transitions on surprisingly short timescales (≲2 months in the rest frame). We find that CL events are preferentially found in lower-Eddington-ratio (f_Edd) systems: Our CL-AGNs have af_Edddistribution that significantly differs from that of a carefully constructed, redshift- and luminosity-matched control sample (Anderson–Darling test yieldingp_AD≈ 6 × 10⁻⁵; medianf_Edd≈ 0.025 versus 0.043). This preference for lowf_Eddstrengthens previous findings of higher CL-AGN incidence at lowerf_Edd, found in smaller samples. Finally, we show that the broad Mgiiemission line in our CL-AGN sample tends to vary significantly less than the broad Hβemission line. Our large CL-AGN sample demonstrates the advantages and challenges in using multi-epoch spectroscopy from large surveys to study extreme AGN variability and physics. 

Abstract Weak emission-line quasars (WLQs) are a subset of type 1 quasars that exhibit extremely weak Lyα+ Nvλ1240 and/or Civλ1549 emission lines. We investigate the relationship between emission-line properties and accretion rate for a sample of 230 “ordinary” type 1 quasars and 18 WLQs atz< 0.5 and 1.5 <z< 3.5 that have rest-frame ultraviolet and optical spectral measurements. We apply a correction to the Hβ-based black hole mass (M_BH) estimates of these quasars using the strength of the optical Feiiemission. We confirm previous findings that WLQs’M_BHvalues are overestimated by up to an order of magnitude using the traditional broad-emission-line region size–luminosity relation. With thisM_BHcorrection, we find a significant correlation between Hβ-based Eddington luminosity ratios and a combination of the rest-frame Civequivalent width and Civblueshift with respect to the systemic redshift. This correlation holds for both ordinary quasars and WLQs, which suggests that the two-dimensional Civparameter space can serve as an indicator of accretion rate in all type 1 quasars across a wide range of spectral properties.